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Abstract 

Exact solutions of nonlinear evolution equations (NLEEs) play a vital role to reveal the internal mechanism of 
complex physical phenomena. In this work, the exact traveling wave solutions of the Boussinesq equation is 
studied by using the new generalized (G'/G)-expansion method. Abundant traveling wave solutions with arbitrary 
parameters are successfully obtained by this method and the wave solutions are expressed in terms of the 
hyperbolic, trigonometric, and rational functions. It is shown that the new approach of generalized (G'/G)-expansion 
method is a powerful and concise mathematical tool for solving nonlinear partial differential equations in 
mathematical physics and engineering. 
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Introduction 

Large varieties of physical, chemical, and biological phe- 
nomena are governed by nonlinear partial differential 
equations. One of the most exciting advances of nonlinear 
science and theoretical physics has been the development 
of methods to look for exact solutions of nonlinear partial 
differential equations. Exact solutions to nonlinear partial 
differential equations play an important role in nonlinear 
science, especially in nonlinear physical science since they 
can provide much physical information and more insight 
into the physical aspects of the problem and thus lead to 
further applications. Nonlinear wave phenomena of 
dispersion, dissipation, diffusion, reaction and convec- 
tion are very important in nonlinear wave equations. 
In recent years, quite a few methods for obtaining 
explicit traveling and solitary wave solutions of non- 
linear evolution equations have been proposed. A variety 
of powerful methods, such as, the inverse scattering trans- 
form method (Ablowitz and Clarkson 1991), the homo- 
geneous balance method (Fan 2000a), the Exp-function 
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method (He and Wu 2006; Akbar and Ali 2012), the 
modified simple equation method (Jawad et al. 2010; Khan 
et al. 2013), the novel (G'/G) -expansion method (Alam 
et al. 2014; Alam and Akbar 2014), the improved (G'/G)- 
expansion method (Zhang et al. 2010), the (G VG)-expan- 
sion method (Wang et al. 2008; Bekir 2008; Zayed 
2009; Zhang et al. 2008; Akbar et al. 2012), the tanh-function 
method (Wazwaz 2005), the extended tanh-function 
method (Fan 2000b; El-Wakil and Abdou 2007), the 
sine-cosine method (Wazwaz 2004), the modified 
Exp-function method (Usman et al. 2013), the gener- 
alized Riccati equation method (Yan and Zhang 2001), 
the Jacobi elliptic function expansion method (Liu 2005; 
Chen and Wang 2005), the Hirota's bilinear method 
(Wazwaz 2012), the Miura transformation method (Bock 
and Kruskal 1979), the new generalized (G '/G)-expan- 
sion method (Naher and Abdullah 2013; Alam et al. 
2013a; Alam and Akbar 2013a; Alam and Akbar 2013b; 
Alam et al. 2013b), the Cole-Hopf transformation 
method (Salas and Gomez 2010), the Adomian decom- 
position method (Adomain 1994; Wazwaz 2002), the 
ansatz method (Hu 2001a; Hu 2001b), the exp(-0(rj))-ex- 
pansion method (Khan and Akbar 2013), the method of 
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bifurcation of planar dynamical systems (Li and Liu 2000; 
Liu and Qian 2001), and so on. 

The objective of this article is to apply the new gener- 
alized (G7G) expansion method to construct the exact 
traveling wave solutions of the Boussinesq equation. 

The outline of this paper is organized as follows: In 
Section Description of the new generalized (GVG)-ex- 
pansion method, we give the description of the new 
generalized (G'lG) expansion method. In Section Ap- 
plication of the method, we apply this method to the 
Boussinesq equation, results and discussions and graphical 
representation of solutions. Conclusions are given in the 
last section. 

Description of the new generalized (G'/G)-expansion 
method 

Let us consider a general nonlinear PDE in the form 



AGG -BGG -EG -C G 



<t>(v, Vt,V x ,V xx ,V t t,V tx , ...). 



(1) 



where v=v(x,t) is an unknown function, <t> is a polyno- 
mial in v(x,t) and its derivatives in which highest order 
derivatives and nonlinear terms are involved and the 
subscripts stand for the partial derivatives. 

Step 1: We combine the real variables x and f by a 
complex variable r] 



v(x, i) = v(rj), rj = x±Vt, 



(2) 



where V is the speed of the traveling wave. The traveling 
wave transformation (2) converts Eq. (1) into an ordin- 
ary differential equation (ODE) for v=v(rj): 



[v,v,v,v", 



(3) 



where i// is a polynomial of v and its derivatives and the 
superscripts indicate the ordinary derivatives with respect 
to n. 

Step 2: According to possibility, Eq. (3) can be inte- 
grated term by term one or more times, yields constant 
(s) of integration. The integral constant may be zero, for 
simplicity. 

Step 3. Suppose the traveling wave solution of Eq. (3) 
can be expressed as follows: 

N N 

vfo) = V a,(rf + M)' + V B t (d + M)-\ (4) 



i=0 



where either a N or B N may be zero, but could be zero 
simultaneously, a, [i=Q,\,2...,N) and /3, (i=l,2,...,N) and d 
are arbitrary constants to be determined and M{rj) is 



M{r,) = (G'/G) 



(5) 



where G=G(rj) satisfies the following auxiliary nonlinear 
ordinary differential equation: 



0. 



(6) 



where the prime stands for derivative with respect to rj; 
A, B, C and E are real parameters. 

Step 4: To determine the positive integer N, taking the 
homogeneous balance between the highest order non- 
linear terms and the derivatives of the highest order 
appearing in Eq. (3). 

Step 5: Substitute Eq. (4) and Eq. (6) including Eq. (5) 
into Eq. (3) with the value of N obtained in Step 4, we 
obtain polynomials in (d+M) N (AT=0,1,2,...) and (d+M)' N 
(A/=0,l,2,...). Subsequently, we collect each coefficient of 
the resulted polynomials to zero, yields a set of algebraic 
equations for a 2 (i'=0,l,2,...,AO and /?,■ {i=l,2,...,N), d 
and V. 

Step 6: Suppose that the value of the constants a, 
(i=0,l,2,...,AO, (i=l,2,...,N), d and V can be found by 
solving the algebraic equations obtained in Step 5. Since 
the general solutions of Eq. (6) are known to us, insert- 
ing the values of a, (z'=0,l,2,...,A/), R t [i=\,2,...,N), d and 
V into Eq. (4), we obtain more general type and new 
exact traveling wave solutions of the nonlinear partial 
differential Equation (1). 

Step 7: Using the general solution of Eq. (6), we have 
the following solutions of Eq. (5): 

Family 1: When B*0, w=A-C and D. -B 2 + 4£(A-C) > 0, 



M(t]) = 



v^Ci sinh ( 


:«') 


+ Ci cosh j 




2cj Cicosh 


(*»: 


) + Ci sinh I 


aw 



(7) 



Family 2: When B*0, co=A-C and O. =BT + 4£{A-C) < 0, 



M( n ) 



B 
2oi 



_q-Ci sin 



J -a . 

2A 1 



■ C 2 cos(^p^ 



2(0 



dcos^pv) +C 2 sin(^v 



(8) 



Family 3: When B * 0, w = A-C and D. = B 2 + 4E{A-Q = 0, 



M{rj) 



G\ B C 2 



(9) 



G I 2co Ci + Civ] 
Family 4: When 5 = 0, co =A-C and A=uE>0, 

MM = (-) = VA ^isinhff >;)+C 2 cosh(f 
' \ G ) w Ci cosh Ifn) +C 2 sinh (^ri) 

(10) 

Family 5: When B = 0, a = A-C and A = w£ < 0, 
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_ A -Ci sin^^J + C 2 cos 
J C lC os(^^) +C 2 sin(^^ 



(11) 



d 4 ^ 2 + 2Bd 3 co + E 2 + B 2 d 2 - 2BdE), n 4 = -(- 8EB 2 o + VV - 
2 V 2 A 4 - 16£V + A 4 -B 4 ),u = A-C,V, d, A, B, C, E are free 
parameters. 

Set 2: K = - ji, (-8£B 2 6) + K 4 A 4 -2K 2 A 4 -16£ 2 <j 2 + A 4 -S 4 ) , 



Application of the method 

In this section, we will put forth the new generalized 
(G7G) expansion method to construct many new and 
more general traveling wave solutions of the Boussinesq 
equation. Let us consider the Boussinesq equation, 



v tt~ v xx~ v xxxx~^{y 2 ^ — 0. 



(12) 



Now, we will use the traveling wave transformation 
Eq. (2) into the Eq. (12), which yields: 



(y 2 -l) v- v m -3{f) 



0. 



(13) 



Eq. (13) is integrable, therefore, integrating with respect 
to T) once yields: 



(V 2 -l) v-v"-3v 2 + K=0, 



(14) 



where K is an integration constant which is to be 
determined. 

Taking the homogeneous balance between highest 
order nonlinear term v and linear term of the highest 
order v" in Eq. (14), we obtain N=2. Therefore, the solu- 
tion of Eq. (14) is of the form: 



v(rj) = a 0 + ai{d + M) + a 2 (d + M) 2 
+ /] 2 (d + M)- 2 , 



-^{d + My 1 
(15) 



where a 0 , a 1( oc 2 , pi, /? 2 and d are constants to be 
determined. 

Substituting Eq. (15) together with Eqs. (5) and (6) 
into Eq. (14), the left-hand side is converted into polyno- 
mials in (d+M) N (N=0,l,2, ) and {d+M)' N (7V=1,2,...). 

We collect each coefficient of these resulted polynomials 
to zero yields a set of simultaneous algebraic equations 
(for simplicity, the equations are not presented) for a 0 , 
a lt a 2t /Si, fi 2 d, K and V. Solving these algebraic 
equations with the help of computer algebra, we obtain 
following: 

Set 1: 



l 

6A 2 



a 0 = — {-A 2 + V 2 A 2 -\2d 2 a 2 + 8Ea-l2Bda-B 2 ) , V=V, 



(17) 



ai =^(2dco 2 +B(j),a 2 = - lj §r,h = 0,/3 2 = 0. Where 
cj = A-C, V, d, A, B, C, E are free parameters. 



Set 3: K 



% «i=0, 



"6 



2u z n 

a2 = -J?^ = Q ^ 2 -*A 2 ^ 



(18) 



where n 5 = {{V 2 - 1)A 2 + 8Eco + 2B 2 ), n 6 = - (16£ 2 w 2 + 
8EB 2 u + B 4 ), n 7 = ((V 2 - 1) 2 A 4 - 256£ 2 w 2 - 128B 2 Eo - 
16B 4 ), o) = A - C,V, A, B, C,E are free parameters. 

For set 1, substituting Eq. (16) into Eq. (15), along with 
Eq. (7) and simplifying, yields following traveling wave so- 
lutions, if Q = 0 but C 2 * 0; C 2 = 0 but Q * 0 respectively: 



1 

6A 2 



B 












cothj 




2a + 


2a 






B 


Vq. 


coth( 


Vn 








2a 


2a 







Hi + 12«2 ( d 



B 






Wet 






tanh( 




2a 


2a 




y ~2A 


B 






Vo 






tanhf 




2a + 


2a 




2A' 



Substituting Eq. (16) into Eq. (15), along with Eq. (8) 
and simplifying, our exact solutions become, if Ci = 0 
but C 2 * 0; C 2 = 0 but C\ *■ 0 respectively: 



1 ,„ / J B V-O fV-O. 

nM = ^ »i + 12« 2 ( d + - + — cot( — n 



2A 



-12k, d 



2a 2a 



a. / v-n 

cot ( rj 



2A 



K 



«4 



12A 4 
2m 



«i 



6A 2 
V = V,d = d. 



, «i = 0,a 2 = 0,/?! 



2« 2 
A 2 "' 



(16) 



where n 1 = {-A 2 + V 2 A 2 - 12d 2 co 2 + 8Eco - \2Bda - B 2 ), 
n 2 = (- 2Edco + 3Bd 2 o + 2c?G) 2 -EB + tfd), n 3 = -(- 2Ed 2 co + 



1 / , B V-Q 

v u (r,)=— 2 U + 12n 2 [d + --—t,n 



2a 2a 



-a 



2A 



B V-Q- fV-O. 
- 12«3 ( d + — tan ( - - rj 



2A 



Substituting Eq. (16) into Eq. (15), together with 
Eq. (9) and simplifying, our obtained solution becomes: 
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Table 1 Comparison between Neyrame et al. (2010) solutions and our solutions 



Neyrame et al. (2010) solutions 



Obtained solutions 



i. If Q = 0 and u(@ = 4v 2 iiq), Case 1 becomes 
v 2l (<D) = -2(X 2 -4/4) coth ; 



* "VC \ 3A 



+ — + a 0 . 



, if A = 1, C= 0, O = X 2 -4/j, 8= 1,f=1,l/=1,f-f = a 0 then the solution is 



v 2l {d>) = -2(A 2 -4f<) coth 2 



/\ 2 -4fiA 3X 



. If Ci = 0 and u(£) = 4i/ 23 (/7), Case 2 becomes: 



v 23 (0) = -2(4M 2 )cot 2 (L ] + - + «„. 



. If = 41/25(1), Case 3 becomes: 



+ f + a 0 - 



|f/\ = 1, C = 0, 0 = A 2 -4fj, 8 = 1, £= 1, V= 1,f -f = 00 then the solution is 



v 23 (0) = -2(4^-A 2 ) cot 2 



/4^-A 2 c 



3A 

+ — + a 0 . 



. If A = 1 , C = 0, S= 1, f = 1,l/= 1,f-f = a 0 then the solution is 



M*)=-2(c^y 



+ f + a 0 . 



= M ni + Un2 ( d+ ^ + c^+c^i) 1 *.<*) = ^ ( * + ( ^ + v coth ( ^ ' 



+12« 3 ( rf + 



£ , C 2 



2a> Ci + C277 



1 I / , VA , /VA 1 ' " 
vi 7 [rj) = —2 < «i + 12« 2 a H tanh — tj 



Substituting Eq. (16) into Eq. (15), along with Eq. 6A 
(10) and simplifying, we obtain following traveling / 









coth 


<a 










'v/A 




coth| 


















tanh 











-1 



^ 2 



wave solutions, if d = 0 but C 2 * 0; C 2 = 0 but d * 0 +12« 3 I d + — tanh I — r] 

respectively: 



A 




Figure 1 Modulus plot singular soliton solution, shape of i/i 6 (n) when A = 4, 8 = 0, C = 1, £= 1, 1/ = 1, d = 0 and -10 < x, f < 10. 
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Substituting Eq. (16) into Eq. (15), together with 
Eq. (11) and simplifying, our obtained exact solutions 
become, if Q = 0 but C2 * 0; C2 = 0 but Q 0 respectively: 



v l8 (7) = ~ «r + 12« 2 ( d + cot ( 7 



6,4 2 



«1 + 12« 2 ^ 



+12« 3 





cot 


V A 




tan 










tan| 




0) 





= -^U 2 (V 2 -l) +2(4fe+£ 2 )-3Dtanh 2 ^^ 



Substituting Eq. (17) into Eq. (15), along with Eq. (8) 
and simplifying yields exact solutions, if C\ = 0 but C2 *■ 0; 
C 2 = 0 but Ci * 0 respectively: 



v 2 M = ^ U^V 2 -1) + 2(4£co + B 2 ) +3Q cot 2 f^y, 



ViAn) = ^|a 2 (V 2 -1) + 2(4fe + J B 2 )-3ntan 2 ('^, 



where n = #-V£. 

Again for set 2, substituting Eq. (17) into Eq. (15), 
along with Eq. (7) and simplifying, our traveling wave 

solutions become, if Ci = 0 but C 2 * 0; C 2 = 0 but Ci * 0 Substituting Eq. (17) into Eq. (15), along with Eq. (9) 



respectively: 



and simplifying, our obtained solution becomes: 
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6A 



\\a 2 (V 1 -1) +2(B 2 + 4£io)-12(o ? C 



Cl + ClY] 



Substituting Eq. (17) into Eq. (15), together with Eq. 
(10) and simplifying, yields following traveling wave solu- 
tions, if Ci = 0 but C 2 * 0; C 2 = 0 but Q * 0 respectively: 



v 2 6 {fj) 



JL ^A 2 (v 2 -l) + (-B 2 + 8Eu) + 12\/A 
x ( B coth( ^ rA -VA coth 2 yj 



V2, 



(^) = ^ ( ^ 2 (^ 2 -l) + {-B 2 + 8£gj) + 12v/A 



x iB cot 



A 



v> 9 (ri)=^(A 2 (V 2 -l) 



x «5tan 



+ VA cot 2 



-fi 2 + 8£w)-12VA 



-v^tan 2 



V 2 7 {l) 



6A Z 



A 2 (V 2 -1) + (-B 2 + 8EcS) + 12y/A 



: | 5tanh( — -rj )-v^tanh 2 ( — ^rj 



where n = x-Vt. 

Similarly, for set 3, substituting Eq. (18) into Eq. (15), 
together with Eq. (7) and simplifying, yields following 
traveling wave solutions, if Q = 0 but C 2 * 0; C 2 = 0 but 
Ci *■ 0 respectively: 



Substituting Eq. (17) into Eq. (15), along with Eq. (11) 
and simplifying, our exact solutions become, if Ci = 0 
but C 2 *■ 0; C 2 = 0 but Q * 0 respectively: 



.^tanh 2 ^ 

n \2A 1 
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V3 2 W=^{« 5 -3Qtanh 2 (-^ 



^coth 2 ^ 

n V 2A ' 



Substituting Eq. (18) into Eq. (15), along with Eq. (8) 
and simplifying, we obtain following solutions, if C l = 0 
but C 2 *■ 0; C 2 = 0 but Q * 0 respectively: 



6A 1 



. , -B a/A , iVA 
h 5 -12w 2 | - — I coth) — /7 



2w cj 



+ 



8co 2 \ 2(0 co I 7 



if o/v-n \ 3« 6 ,/V-n 



« 5 + 3D tan 2 ' - > 



2A n ) Q. COt I 2A n 



, -5 a/A , A/A 
M 5 -12cy | — + — tanhj — ^ 



2w cj 



/I 



8a 2 2a 6) \ A 



Substituting Eq. (18) into Eq. (15), along with Eq. (11) 
Substituting Eq. (18) into Eq. (15 , along with Eq. (9) , . ,., . , t . , / , J? i_ 'c 

, . ... . , . , ,. , and simplifying, our obtained exact solutions become, it 



and simplifying, our obtained solution becomes: 



Ci = 0 but C 2 * 0; C 2 = 0 but Q * 0 respectively: 



Substituting Eq. (18) into Eq. (15), along with Eq. (10) 
and simplifying, yields following exact traveling wave solu- 
tions, if Ci = 0 but C 2 *■ 0; C 2 = 0 but Q * 0 respectively: 



V3 8 (7)= 6 ^^ 5 -12 W ■ . o 



- + - 



cot 



3« 6 / -B V-A / \/-A x 

+8^ 2w + ^r C0t h^' 7 , 
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V3,{l) = 




where rj = x-Vt. 
Results and discussions 

It is worth declaring that some of our obtained solu- 
tions are in good agreement with already published 
results which are presented in the following tables 
(Table 1). 

Beside this table, we obtain further new exact traveling 
wave solutions v 2l (rj), v 2l {rj), v 26 (rj)-v 29 (ij), v^itf-v^rj), 
V3 1 (ij)-V3 9 (ij) in this article, which have not been reported 
in the previous literature. 

Graphical representation of the solutions 

The graphical illustrations of the solutions are given 
below in the figures with the aid of Maple (Figures 1, 2, 
3, 4 and 5). 

The solutions corresponding to vi^ij), V\ 3 (rf), Vi 5 (tj), 

Vu(l)l V 2Al)> V 2s(l)i V 31 (V), v 3M- V 3Al) 

is identical to the solution Vi 6 (tj) , the solution corre- 
sponding to Vi 2 (rj) is identical to the solution v 27 (rj), the 
solution corresponding to v 2l (rf) is identical to the solu- 
tion Vi 7 (rj), the solution corresponding to v 2s {rf), v^rj), 
V3i(tj), v 3g(t]) v 39 (tj) is identical to the solution v 3a (t]) 



and the solution corresponding to Vi^rj), Vi s {rj), Vi 9 (fj), 
v 23 (rj), V 2i {rf) is identical to the solution v 2ci (rj). 

Conclusion 

In this paper, we obtain the traveling wave solutions 
of the Boussinesq equation by using the new ap- 
proach of generalized (G'/G) -expansion method. We 
apply the new approach of generalized (G'/G)-expan- 
sion method for the exact solution of this equation 
and constructed some new solutions which are not 
found in the previous literature. This study shows 
that the new generalized (G'/G)-expansion method is 
quite efficient and practically well suited to be used 
in finding exact solutions of NLEEs. Also, we observe that 
the new generalized (G7G)-expansion method is straight- 
forward and can be applied to many other nonlinear 
evolution equations. 
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